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ABSTRACT 

We investigate the viability of the braneworld model with energy exchange be- 
tween the brane and bulk, by using the most recent observational data related to the 
background evolution. We show that this energy exchange behaves like a source of dark 
energy and can alter the profile of the cosmic expansion. The new Supernova Type 
la (SNIa) Gold sample, Supernova Legacy Survey (SNLS) data, the position of the 
acoustic peak at the last scattering surface from the Wilkinson Microwave Anisotropy 
Probe (WMAP) observations and the baryon acoustic oscillation peak found in the 
Sloan Digital Sky Survey (SDSS) are used to constrain the free parameters of this 
model. To infer its consistency with the age of the Universe, we compare the age 
of old cosmological objects with what computed using the best fit values for the 
model parameters. At 68% level of confidence, the combination of Gold sample SNIa, 
Cosmic Microwave Background (CMB) shift parameter and SDSS databases provide 
Q m = 0-29±g:gl, n A = -0.71 and [i = -0.40ig|{j, hence a spatially flat Universe 
with n K = o.oot°;ot The same combination with SNLS supernova observation give 
fi m = 0.27t°;o2, M A = -0.74t°;p2 and /i = O.OOt^o consequently provides a spatially 
flat Universe fl K = -0.01±g;^. These results obviously seem to be compatible with 
the most recent WMAP results indicating a flat Universe. 

Key words: Cosmology - methods: numerical - methods: statistical - cosmology: 
theory - cosmology: cosmological parameters - cosmology: early Universe - cosmology: 
observations 



1 INTRODUCTION 

Recent observations of type la supernovas (SNIa 
sugge st the expansio n of the Universe is accelerat 



ing (lA.G. Riess el all Il998l; IS~Perlmutter et all Il999l; 
lA. G Riess el. all 12004 |j. L. Tonrv el. all I2003I V As it 
is well known all usual types of matter with positive 
pressure generate attractive forces, which decelerate the 
expansion of the Universe. A "dark energy" component 
with negative pressure was suggested to account for the 
invisible fuel that drives the current acceleration of the 
Universe. Although the nature of such dark energy is 
still speculative, an overwhelming flood of papers has 
appeared which attemp t to describe it by devising a great 
varie t y of models (see (|V. Sahni el. al.l l2000l: IS. Weinberg 
1 19891 : I J. A. S. Lima! 12004 |E. J. Copeland el. al.l l200rj ; 
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IC. Armendariz-Picon elaU |2000| : |j. S. Bagla el. al.l 120031 ) 

for recent reviews) . Available models of dark energy differ in 

the value and variation of the equation of state parameter, 

w, during the evolution of the Universe. Among them are 

cosmological constant A, an evolving scalar field (referred 

to by some as quintessence), the phantom energy, in which 

the sum of the pressure and energy density is negative , the 

quintom model, the holographic dark energy , the Chap- 

lygin gas, and the Cardassion model. Another approach 

dealing with this problem is using the modified gravity 

by changing the Einstein-Hilbert action. Some of models 

as 1/72 and logarithmic models provi de an acceleration 

for the Universe at the present time (IS. Weinberg 19891 : 
„ , , ^ rnu™-,i 1^ t „ — r; — 1™ , , , — 1 



M.Ca rroll 200 il:lP. J. 



20031: IC. L. Bennett eLa! 



E. Peebles 2003; T. Padmanabhan 
20031: IH.V. Peiris el. al.l 12003 



D. N Spergel el. al.l l2003al; I L F Miranda et all I2OO1I 
S. Rahvar el. all 120071: IC. Wetterichl ll99Sl; Ib. Ratra et. al 



19881 : I J. A. Frieman et. al.l Il995l : I M S. Turner et. al 
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19971; I R. R. Caldv u l] . i . a] 
I. Zlatevl Il999l: 



2002 



2003 



2000 



19981 ; I A. R. Liddld Il998l: 



P. J. 



IP. J. E. Peebles el 
S. Arbabi-Bidgoli cl. al 



Stei nhardt 



1988 



19991: 



D. F. Torres 



M. 



■S. Perlmutter el al.l Il999t . 
Doran el. all l200ll; |R. R. Caldwell el. al 
200: 



.R R. Caldwell el. al 
20061; lL. Wang el. al 
|L. Page el. al l 12003 



R. R. Cal dwell 
IT 



2004 
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2003 



2003 
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Franca et. al.l 
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2006a; B. Wang et. al 
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Independent of these challenges, we deal with the dark 
energy puzzle. In recent years, theories of large extra dimen- 
sions, in which the observed Universe is realized as a brane 
embedded in a higher dimensional spacetime, have received 
a lot of interest. According to the braneworld scenario, the 
standard model of particle fields are confined to the brane 
while, in contras t, the gravity is free to propagate in the 
whol e spacetime ||L. Randall el. ailll999l ; |G. R. Dvali et. ail 
2000). In these theories the cosmological evolution on the 
brane is described by an effective Friedmann equation 
that incorporat es non-trivially with the effects of the bulk 
into th e brane d P. Binetruv el."aT~ll2000l; |A Sheykhi el. all 
l2007al ; |A. Sheykhi. B. Wang and R.G. Cail l2007d ). An 



interesting consequence of the braneworld scenario is that 
it allows the presence of five-dimensional matter which 
can propagate in the bulk space and may interact with 
the matter content in the braneworld. It has been shown 
that such interaction can alter the profile of the cosmic 
expansion and leads to a behavior that would resemble the 
dark energy. The cosmic evolution of the braneworld mod- 
els with energy exchange between the br ane and bulk has 
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In the framework of the braneworld scenarios, many at- 
tempts to observationally detect or distinguish brane effects, 
on the evolution of our Universe, from the us ual dark energy 
phys i cs have been discussed in th e lite rature (IN. Pires el. all 
20061; IS. Capozziello et. aill2004l ). In (|V. Sahni el. al.ll2002l ; 



V. Sahni et. aill2003l ) a class of braneworld models has been 



investigated. A new and interesting feature of this class of 
models is that the acceleration of the Universe may be a 
transient phenomenon, which cannot be achieved in the con- 
text of our current standard scenario, i.e., the ACDM model 
but could reconcile the supernova evidence for an accel- 
erating Universe with th e requirements of string/M-theory 
I W. Fischler el. alJlioOlh . The purpose of the present work 
is to disclose the effect of energy exchange between the brane 
and bulk in Randall- Sundrum II braneworld scenario on the 
evolution of the Universe. Giving the wide range of cosmo- 



logical data available, we are able to test the viability of this 
class of braneworld models by putting recent observational 
constraints on its free parameters. 

We have three independent types of observational 
constraints for the dark ener gy models: (i) the 
pernova distanc e 



A. G. Riess el. ail 



mod ulus ( A.G. Riess el. al.l 

20041 ; IS. Perlmutter el. al. 



1998 



1998 



S. Nobili et. alj 12005), (it) the dynamica l evidence for 



matter density (D^N. Spcrg el el. al 



L. Knox el. al.l 



200381) and (Hi) the 
200ll ; I Hu. W. et. all 



age o f the Universe 

l200ll ). Besides, a great success has been scored in high preci- 
sion measur ements of CMB anisotropy, as well as in galaxy 
clust e ring (ID. H. Weinberg el. ail 120051; IU.Seliak et.~al] 



120041 ; |A. Refregieij|2003l ; |C. Hevmans et. ail 120051 ). Among 
these observations, the age of the Universe is one of the 
most pressing pieces of data disclosing information about 
dark energy. Indeed, any limit on the age of the Universe 
during its evolution with redshift will reveal the nature 
of dark energy. This is due to the fact that dark energy 
influences the evolution of the Universe. However, different 
models of dark energy may lead to the same age of Universe 
at z = 0. To lift this degeneracy, we should examine the 
age of the Universe at different stages of its evolution and 
compare it with the estimated age of high-redshift objects. 
This procedure constrains the age at different stages, being 
a powerful tool to test the viability of di fferent models 
||B. Wang et. alll2006bl ; |A. Friaca el. al.ll2005l ). 

This paper is organized as follows: In section 2, we in- 
troduce a braneworld model with energy exchange between 
the brane and bulk, the cosmology of this model, its free 
parameters and background dynamics of the Universe gov- 
erned by the effective Friedmann equation. We also show 
how this model can exhibit acceleration expansion of our 
Universe. Most limitations regarded to this interaction in 
our model are introduced. We investigate the geometrical 
effects of underlying braneworld cosmology in section 3. In 
section 4, we test the viability of our model by putting some 
constraints on the parameters of the model. For this aim, we 
use the new Gold sample and Legacy Survey of SNIa data 
(Riess et al. 2004; Astier et al. 2005), its combination with 
the position of the observed acoustic angular scale on CMB 
and the baryonic oscillation length scale. In section 5, we 
compare the age of the Universe in this model with the age 
of old cosmological objects. The last section is devoted to 
conclusions and discussions. 



2 BRANEWORLD WITH BRANE-BULK 
ENERGY EXCHANGE 



We start from the following action 
S = i / dS W=ff(-R-2A) + 



j- m 



(i) 



where R is the 5D scalar curvature and A < is the bulk 
cosmological constant, g and g are the bulk and the brane 
metrics, respectively. Throughout this paper we choose the 
unit k 2 = 1 as the gravitational constant in five dimension. 
We have also included arbitrary matter content both in the 
bulk and on the brane through L™ ulk and L™ rane respectively. 
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a is the positive brane tension. The field equations can be 
obtained by varying the action, equation ([1]), with respect 
to the bulk metric qab- The result is 



Gab + A.gAB = Tab- 



(2) 



For convenience we choose the extra-dimensional coordinate 
y such that the brane is located at y = and bulk has Z2 
symmetry. We are interested in the cosmological solution 
with a metric 

ds 2 = —n 2 (t,y)dt 2 +a 2 (t,y) , yijdx l dx : '+b 2 (t,y)dy 2 , (3) 

where 7^ is a maximally symmetric 3-dimensional metric for 
the surface (£=const., y=const.), whose spatial curvature is 
parameterized by K = — 1, 0, 1. The metric coefficients n and 
b axe chosen n(t, 0) = 1 and b(t, 0) = 1, where t is cosmic 
time on the brane. The total energy-momentum tensor has 
bulk and brane components and can be written as 



Tab = Tal 



+T Ai 



+T Ai 



(4) 



The first and the second terms are the contribution from the 
energy-momentum tensor of the matter field confined to the 
brane and the brane tension 



T A B \ b rane = diag(-/0, p, p, p, 0) — — , 



rrtA 

1 B 



diag(- 



-(7,0) 



(5) 
(6) 



where p, and p, being the energy density and pressure on 
the brane, respectively. In addition, we assume an energy- 
momentum tensor for the bulk content with the following 
form 



1 o 



B \bulk — 








(7) 



The quantities which are of interest here are T% and T° 5 , 
as these two enter the cosmological equations of motion. In 
fact, T° 5 is the term responsible for energy exchange between 
the brane and the bulk. Integrating the (00) and the (if) 
components of the field equations ((2)1 across the brane and 
imposing Z2 symmetry, we have the jump across the brane 

£ = -~(P + <7). (8) 
ao b 

^ = i(2p + 3p-a), (9) 
no b 

where 2a' + = —2a'_ and 2n' + — —2n'_ are the discontinuities 
of the first derivative and primes denote derivatives with 
respect to y. In addition, as usual, the subscript " 0" denotes 
quantities are evaluated at y = 0. 

Substituting the junction conditions i.e. equations ((8} 
and © into the (55) and (05) components of the field equa- 
tion ([5]), we obtain the modified Friedmann equation and 
the semi-conservation law on the brane 
1 

~36 



2H 2 + H+* = -^-[a(3p-p)+p(p + 3p)] 
"0 



p + 3H(p + p) = -2T° 5 , (11) 

where H — do /ao is the Hubble parameter on the brane and 
dots denote time derivative. We shall assume an equation 



of state p = wp which represents a relation between the 
energy density and pressure of the matter on the brane. 
The bulk matter contributes to the energy content of the 
brane through the bulk pressure terms T% and T%. In order 
to derive a solution that is largely independent of the bulk 
dynamics, we should neglect T 5 term by assuming that the 
bulk matter relative to the bulk vacuum energy is much less 
than th e ratio of the brane m atter to the brane vacuum 
energy (|E. Kiritsis el. al"1l2003T ). Considering this we get 



2H 2 + H + ^ = IP (1 - 3w) - (3p 2 (1 + 3w) + ^, (12) 



p + 3Hp(l + w) = -2T° 5 , 



(13) 

where we have used the usual definition /3 = 1/36, 7 = Pa 
and A = (A + a 2 /6). Assuming the Randall- Sundrum fine- 
tuning A = A + a 2 /b — holds on the brane, one can easily 
check that the Friedmann equation (I12|l is equivalent to the 
following equations 

H 2 + ^ = f3p 2 + 2 1P + x, 



X + 4H X = 4T 5 (f3p + 1 ). 



(14) 
(15) 

Equation (|14|l is the modified Friedmann equation describ- 
ing cosmological evolution on the brane. The auxiliary 
field x incorporates non-trivial contributions of dark energy 
which differ from the standard matter fields confined to the 
brane. It is worth noting that the flow of the mass-energy 
from the bulk onto the brane may resemble as the dark en- 
ergy. Indeed it can influence the backgroun d evolution of the 
Univ erse and leads to acceleration (see e.g. (|E. Kiritsis el. all 
2003)). One may argue that whether the energy exchange 
between the brane and bulk becomes dark matter or not? 
To answer to this question, one should consider an inter- 
action between dark matter and dark energy on the brane 
which is not clear yet. Besides in order to have the equation 
of state in the bulk, a particular model of the bulk matter 
is required which is not cl ear yet, because we do no t exactly 
know the bulk geometry (|C. Bogdanos el. al.ll2007l ). So now 
in our coarse-grained model we ignored this effect. 

We are also interested in the scenarios where the energy 
density of the brane is much lower than the brane tension, 
namely p a, therefore equations (|14|l and (|15p can be 
simplified in the following form 



rr2 K 

at 



27P + X, 



(16) 







X + AH X = 4 7 T° 5 . (17) 

Then we take ansatz T% = AHa? f or the brane-bulk en- 
ergy exchange (|R.G. Cai el. al.ll2006h . where A and p, are 
arbitrary constants and thereafter we have omitted the "0" 
subscript from the scale factor on the brane for simplicity. 
For this ansatz, one can easily check that equation (|17|l has 
the following solution 



C , 4 7 ^1 
— r H or 

a 4 p + 4 



(18) 



where C is an integration constant usually referred to the 
dark radiation term. In a similar way, inserting T° 5 into 
equation (|13|) . we get 

2A „ 



P - 



P<> 



/' - 



(19) 
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where po is the present matter density of the Universe with 
equation of state w — 0. Finally, inserting p and x i n t° 
equation (I16|) . we obtain the modified Friedmann equation 
on the brane 



H 1 



9,-kG? 



2A 



(p + 3)(/x + 4) 



K 
7a 



(20) 



where Gn = S-y/A-K is the 4D Newtonian constant, p m — 
poa~ 3 is matter energy density and we have neglected the 
dark radiation term ~ a -4 , namely C — 0, because we are 
more interested in the prob of late time era. Using the value 
of present critical density, 

P° = ^Z7^' ( 21 ) 

the effective Friedmann equation in terms of dimensionless 
quantities and redshift parameter 1 + z = a~ can be written 

as 

H 2 = H%[n m (l + z) 3 -n A {l + z)->*-({l tot -l)(l + z) 2 ],(22) 



where 

Qm — 

Q A = 

"tot = 



PO 



0/ 



2.4 



3K 



StvGnPc 



p c (p + 3)(p + 4)' 
Q m — Ha = 1 + Ok- 



(23) 

(24) 
(25) 



As one can see from equation (|22[) . the free parameters of 
this model are very similar to those of ACDM models, but 
this is quite accidental and is due to our specific ansatz for 
the energy exchange term T° 5 . Indeed, the energy exchange 
term Qa in equation (|22l) which behaves such as cosmolog- 
ical constant term in ACDM models is originated from the 
bulk matter content (see equation Q). This is completely 
different from the origin of the corresponding term in the 
ACDM models. In order to get the late time acceleration ex- 
pansion profile for the Universe, this term plays a crucial role 
here. Therefore the braneworld model with energy exchange 
between the brane and bulk, gives a very useful framework 
for comparing the ACDM general relativistic cosmology to 
a modified gravity alternative. 

To see how our model can exhibit acceleration expan- 
sion of our Universe, we study the behavior of the acceler- 
ation parameter. One can easily show that the acceleration 
parameter in this model can be written as 

q ee —-^ = -| [Q m (l + z) z + Q A {p + 2)(1 + 2)-"] . (26) 

As it can be seen in figure [1] increasing p causes the Uni- 
verse to accelerate earlier. In figure[2]we compare this model 
and ACDM just according to the acceleration parameter. In- 
creasing (decreasing) the value of Ha (Q\ in ACDM model) 
causes the Universe to enter accelerating epoch earlier. As 
we will see in the following section using the best fit values 
for model parameters, acceleration parameter in the present 
time at la confidence level is q(z = 0) = 0.42t(j'i9 while for 
ACDM model is q(z = 0) = 0.59 j:^. 

Now an interesting question that arises is: can this 
model predict dynamics of the Universe? In other words, 
For what values of the free parameters, theoretical model is 
consistent with the observational tests? 

In the forthcoming sections we will see what constraints 
to the described model are set by recent observations. As 




Figure 1. Acceleration parameter (g = d/a_ffg) in the 
braneworld model as a function of redshift for various values of 
fjj. Here we chose = 0.0, f2 m = 0.30 and Qa — —0.70. 



a matter of fact we examine the free parameters of model 
more carefully. Indeed we let the parameters scan their phase 
space and using likelihood statistics, the best fit values which 
maximize likelihood function will be retrieved. 



3 GEOMETRICAL EFFECTS OF 
BRANEWORLD MODEL 

The cosmological observations are mainly affected by the 
background dynamics of the Universe. So before starting 
some main observational tests to explore braneworld cosmol- 
ogy we investigate how the free parameters of this model al- 
ter the background dynamics by using the measurable quan- 
tities introduced in this section. We believe they give deep 
insight throughout this model. For this purpose, we study 
the effect of the braneworld model on the geometrical pa- 
rameters of the Universe all together. 



3.1 comoving distance 

The radial comoving distance is one of the basic parameters 
of cosmology. For an object with the redshift of z, using the 
null geodesies in the FRW metric, the comoving distance is 
obtained as: 



r(z;U m ,Q, A ,fJ>) 



H y/\UV\ 



dz' 



H(z')/H 



(27) 



where 



(x, sin a:, sinh x) for K = (0, 1, — 1) 



(28) 



and H(z; Sl m , ^a, p) is given by equation (|22|l . 

By numerical integration of equation (|27[) . the comoving 
distance in terms of redshift for different values of p is shown 
in figure [3] Increasing the p results in a longer comoving 
distance. According to this behavior by fine tuning the value 
of p in addition to Qa and f2 m , one can expect to explain 
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Figure 2. Upper panel shows acceleration parameter (q = 
a/H^a) in the braneworld model as a function of scale factor for 
various values of Qa an d A 1 = —0.40. Lower panel corresponds to 
the same function for the flat ACDM. We chose the fiat Universe. 



Figure 3. Comoving distance, r(z; Q m , fi) (in unit of c/Hq) 
as a function of redshift for various values of /i. We fixed fi^ = 
0.0, Q m = 0.3 and tt A = -0.7. 




the observational results given by supernova as a standard 
candle to measure distance in the observational cosmology. 

3.2 Angular Size 

The apparent angular size of an object located at the cos- 
mological distance is another important parameter that can 
be affected by the cosmological model during the history of 
the Universe. An object with the physical size of D is related 
to the apparent angular size of 9 by: 



D = d A e, 



(29) 



where cLa = r(z; Q. m , Qa, /x)/(1 + z) is the angular diameter 
distance. The main applications of equation (|29[) is on the 
measurement of the apparent angular size of acoustic peak 
on CMB and baryonic acoustic peak at the high and low 
redshifts, respectively. By measuring the angular size of an 
object in different redshifts (the so-called Alcock-Paczynski 
test) it is possible to pr obe the validity of braneworld model 
|C. Alcock el. all 1 19T9T ) . The variation of apparent angular 
size A8 in terms of Az is given by: 



Figure 4. Alcock-Paczynski test comparing Az/A8 as a function 
of redshift for five different fi normalized to the case with f2 m = 
0.30 and Q A = -0.70 and n = 0.0 (fiat Universe fl K = and 
ACDM ). 



Az 

— = H(z;Qm,QA, n)r(z;Q. m ,SlA, (30) 

Figure [4] shows Az/A6 in terms of redshift, normalized to 
the case with Q m — 0.30, Qa = —0.70 and /i = 0.0 (flat Uni- 
verse, fix = 0.0). The advantage of Alcock-Paczynski test 
is that it is independent of standard candles and a standard 
ruler such as the size of baryonic acoustic peak can be used 
to constrain the braneworld model. 



3.3 Comoving Volume Element 

The comoving volume element is another geometrical pa- 
rameter which is used in number-count tests such as lensed 
quasars, galaxies, or clusters of galaxies. The comoving vol- 
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4.1 Supernova Type la: Gold and SNLS Samples 




Figure 5. The comoving volume element in terms of redshift 
for various fj. exponent. Increasing fi shifts the position of the 
maximum value of volume element to higher redshifts. We fixed 
n K = 0.0. 



Parameter 


Prior 




f2tot = n m — &a 




Free 




[0.00,1.00] 


Top hat 


a b h 2 


0.020 ± 0.005 


Top hat (BBN) 


n A 


[-3.00,1.00] 


Top hat 


h 




Free 


n 




Free 



Table 1 ■ Priors of the parameter space, used in the likelihood 
analysis dA. Melchiorri el. al.|[20o3; IW. L. Freedman el. aLll200lj; 
IX. Zhang el. alj|2005l) . 



ume element in terms of comoving distance and Hubble pa- 
rameter is given by: 



dV 
dzdQ, 

r 2 {z; O m , Oa, (j)/H{z; O m , Oa, /i)(31) 



According to figure[5] the comoving volume element becomes 
large for larger value of fi in the flat Universe. 



The Supernova Type la experiments provided the main 
evidence of the existence of dark energy. Since 1995 two 
teams of the High-Z Supernova Search and the Super- 
nova Cosmology Project have dis covered several type la su- 
pernovas at the high red shifts (IS. Perlmutter el. al.l 1 19991 : 



IB. P. Schmidt el. alj[l99l ') . Recently, Riess et al.(2004) have 



announced the discovery of 16 type la supernova with 
the Hubble Space Telescope. They determined the lumi- 
nosity distance of these supernovas and with the previ- 
ously reported algorithms, obtai ned a uniform 157 Gold 
sample of type la supernovas (lA. G. Riess el. all 120041 : 
I J. L. Tonrv el. al.l 120031 : IB. J. Barris el. al.l 12004 ). Recently 
a new data set of Gold sample with smaller system- 
atic error containing 15 6 Supernova la has been released 
l|The Gold datasetl l2006). In this work we use this data set 
as new Gold sample SNIa. 

More recently, the SNLS collaboration released the 
first ye ar data of its planne d five-year Supernova Legacy 
Survev (|P. Astier el. al.ll2005l ). An important aspect to be 
emphasized on the SNLS data is that they seem to be in a 
bett er agreement with WMA P results than the Gold sam- 
ple (|H. K. Jassal el. al.ll20Q6h . We compare the predictions 
of the braneworld model for apparent magnitude with new 
SNIa Gold sample and SNLS data set. The observations of 
supernova measure essentially the apparent magnitude m 
including reddening, K correction, etc, which are related to 
the (dimensionless) luminosity distance, Dl, of an object at 
redshift z through 



m = M + 5log D L (z; Q m ,O.A,n), 



(32) 



where 

D L {z; O m , Oa, 
Also 



M 



M + 5 log 



(34) 



where M is the absolute magnitude. The distance modulus, 
5ft, is defined as 



5ft = m — M — 5 log DL,(z;Q.m,Q a, 

1 Mpci +25 ' 



+5 log 



(35) 



4 OBSERVATIONAL CONSTRAINTS ON THE 
MODEL USING BACKGROUND 
EVOLUTION OF THE UNIVERSE 

In this section, at the beginning, we examine braneworld 
model by SNIa Gold sample and supernova Legacy Survey 
data. Then to make the model parameter intervals more con- 
fined, we will combine observational results of SNIa distance 
modules with power spectrum of cosmic microwave back- 
ground radiation and baryon acoustic oscillation measured 
by Sloan Digital Sky survey. Table 1 shows different priors 
of the model parameters used in the likelihood analysis. 



5ft = 5\og D L (z;fl m ,n A ,fi) + M. (36) 

In order to compare the theoretical results with the obser- 
vational data, we must compute the distance modulus, as 
given by equation (|35p . For this purpose, the first step is to 
compute the quality of the fitting through the least squared 
fitting quantity x 2 defined by: 

Y 2 (M, 0„,, Q A , A*) = 

[5ftobsQz») -5ft th (zi;O m ,OA,M,^)] 2 ,o 7 x 

i 

where Oi is the observational uncertainty in the distance 
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modulus. To constrain the parameters of model, we use the 
Likelihood statistical analysis: 



(38) 



where N is a normalization factor. The parameter M is a 
nuisance parameter and should be marginalized (integrated 
out) leading to a new \ 2 defined as: 



X 2 = -2m 



c(M,n m ,n A ,n)dM. 



(39) 



Using equations (|37|l and (I39|) . wc find 

x(Q. m ,£lA,y) = x 2 {M = 0, n m , Q. A , y) 

B(n m ,nA,y) 2 



c 



+ ln(C/27r) 



where 



B(Q m , Qa, y) = 

[K obs (zi) - 3? t h(zi;n m ,n A ,/i,M = 0)] 



E 



and 



(40) 



(41) 



(42) 




0.2 0.4 0.6 



Figure 6. Distance modulus of the SNIa new Gold sample in 

terms of redshift. Solid line shows the best fit values with the 

+0.10 



corresponding parameters of Hg = 63.66, fl n 



0.51 



n A 



-0.75±? •?? and fi 



with xliJN, 



d.o.f 



0.76^'gg in lcr level of confidence 
0.92 for brane world model. 



Equivalent to marginalization is the minimization with re- 
spe ct to M. One can show that \ 2 can be expanded in M 
as f|S. Nesseris el. al.ll2004T ) 



2 

XSNIe 



(Q, m ,Q,A,y) = X (M — 0, fi m , &a, A*) 
-2A?B + M 2 C, 



which has a minimum for M = B/C: 

XSNIaG^m, ^A, A*) = (M = 0, fi m , f2 A , ft) 

B(n m ,04, a*) 2 



c 



(43) 



(44) 



Using equation (|44|l we can find the best fit values of model 
parameters as the values that minimize XsNia(^™> Qa, £*)• 
The best fit values for the parameters of model by using 
supernova data are f2„ 



0.51^30, «A 



-0.75l?; 4 2 and 



H = 0.76li95 with XTnin/Nd.o.f = 0.92 at lcr level of confi- 

+ ?;^|. The best 



dence. These values imply that Ok = —0.26 
fit values for the parameters of model by using SNLS su- 
pernova data are f2 m = 0.06^o'q 6 , Qa = — 1.84^ ' S g and 
(j, = -1.34to.?o with X 2 min/Nd.o.f = 0.87 at la level of con- 
fidence. The corresponding value of Qk at la confidence 



level is Q,k 



-0.90 i 



Figures [6] and [7] show the com- 



parison of the theoretical prediction of distance modulus by 
using the best fit values of model parameters and observa- 
tional values from new Gold sample and SNLS supernova, 
respectively. Figures [5] and show relative likelihood for free 
parameters of brane model. 

4.2 Combined analysis: SNIa+CMB+SDSS 

To obtain more confined acceptable intervals of model free 
parameters, now we combine SNIa data (from SNIa new 
Gold sample and SNLS) with CMB data from the WMAP 
and recently observed baryonic peak from the SDSS. We also 
examine the peaks positions of power spectrum in addition 
to the common shift parameter. 




Figure 7. Distance modulus of the SNLS supernova data in 
terms of redshift. Solid line shows the best fit values with the 



corresponding 

n A = -i.84; 

with X 2 mi JN d . .f 



parameters of Hq = 69.38, Q m = 0.06^Q gg, 
1.34j] ' 10 in lcr level of confidence 
0.87 for brane world model. 



0.59 and A 4 



Before last scattering, the photons and baryons are 
tightly coupled by Compton scattering and behave as a fluid. 
The oscillations of this fluid, occurring as a result of the bal- 
ance between the gravitational interactions and the photon 
pressure, lead to the familiar spectrum of peaks and troughs 
in the averaged temperature anisotropy spectrum which we 
measure today. The odd and even peaks correspond to max- 
imum compression of the flu id and to rarefaction, respec- 
tively i| Wayne Hu el. alJI 1997ft . In an idealized model of the 
fluid, there is an a nalytic relation for the location of the m-th 
peak: l m w ml a (|W. Hu el. al.lll995l : fHu. W et. alJl200lft 
where I a is the acoustic scale which may be calculated ana- 
lytically and depends on both pre- and post-recombination 
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Figure 8. Marginalized likelihood functions of three parameters 
of model (O m , VI a and fi). The solid line corresponds to the like- 
lihood function of fitting the model with SNIa data (new Gold 
sample), the dashdot line with the joint SNIa+CMB data and 
dashed line corresponds to SNIa+CMB+SDSS. The intersections 
of the curves with the horizontal solid and dashed lines give the 
bounds with lcr and 2cr level of confidence respectively. 



Figure 9. Marginalized likelihood functions of three parameters 
of a braneworld model (f2 m , !lyi and fi). The solid line corresponds 
to the likelihood function of fitting the model with SNIa data 
(SNLS), the dashdot line with the joint SNIa+CMB data and 
dashed line corresponds to SNIa+CMB+SDSS. The intersections 
of the curves with the horizontal solid and dashed lines give the 
bounds with lcr and 2a level of confidence respectively. 



physics as well as the geometry of the Universe. The acous- 
tic scale corresponds to the Jeans length of photon-baryon 
structures at the l ast scattering surface some ~ 379 Kyr 
after the Big Bang (|D. N. Spergel el. alj|2003af ). The appar- 
ent angular size of acoustic peak can be obtained by dividing 
the comoving size of sound horizon at the decoupling epoch 
r B (zdec) by the comoving distance of observer to the last 
scattering surface r(zdec) 



The size of sound horizon at the numerator of equation (|45|l 
corresponds to the distance that a perturbation of pressure 
can travel from the beginning of the Universe up to the last 
scattering surface and is given by 
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Figure 10. Constant acoustic angular scale in the joint space 
of Q m and (upper panel). Lower panel shows dependence of 
acoustic angular scale on the Q m and cosmological constant. 



where v s (z)~ 2 — 3+9/4x pb{z) / p rad (z) is the sound velocity 
in the unit of speed o f light from the bis bang up to the 
last s cattering surface (|M. Doran el. alJbOQll ; fW. Hu el. al.l 
1 19951 ) an d the redshift of the last scattering surface, Zdec, is 
given by (|W. Hu el. al.lll995h : 



Zd 



!C = 1048 [l + 0.00124(, 



[l+gi(u m ) n ], 



gi = 0.0783(w 6 )- a238 [l+39.5(^) a763 l 1 



.92 



0.560 1 + 21.1(^6 



1.81T - 1 



(47) 



where ui m = flmh 2 , U)b = flbh 2 and p rad is the radiation den- 
sity, fib is relative baryonic density to the critical density at 
the present time. Changing the parameters of the model can 
change the size of apparent acoustic peak and subsequently 
the position of I a = K /9 a in the power spectrum of temper- 
ature fluctuations at the last scattering surface. The simple 
relation l m ~ ml a however does not hold v ery well for the 
peak s although it is better for hig her peaks ll Hu. W. et. al.l 
l200ll : I Michael Doran el. alj|200ll ). Driving effects from the 



decay of the gravitational potential as well as contributions 
from the Doppler shift of the oscillating fluid introduce a 
shift in the spectrum. A good parameteriza tion for the lo- 
cation of the peaks and troughs i s given by l| Hu. W. et. al.l 
l200ll ; I Michael Doran el. al.ll200ll ) 



lm — lA(m — 4>m), 



(48) 



where <j>m is phase shift determined predominantly by pre- 
recombination physics, and are independent of the geom- 
etry of the Universe. The location of acoustic peaks can 
be determined in model by equation l|4"8"|l w ith 4> m (uo m ,u)b)- 
Doran et. al. (| Michael Doran el. al.l [2001). have recently 
shown that the first and third phase shifts are approximately 
model independent. The values of these shift parameters 
have b een reported as: < bi(u)m , ^b) — 0-27 and <j>3(uJm , Ub) 



0.341 (I Hu. W. et. all boOll; Iw. J. Percival et. al. I |2002l ; 



I Michael Doran el. al.ll200ll ). According to the WMAP ob- 
servations: h = 220.1 ± 0.8 and Z3 = 809 ± 7, so the corre- 
sponding observational values of read as: 



7obs I 

l A \h = 



h 



A \h 



(1 -4n) 
h 



= 299.45 ± 2.67, 
304.24 ± 2.63, 



(3 - 0s) 

their Likelihood statistics are as follows: 



2 

Xh 



and 



Xh 



[/obs I ;thi 
'a \h ~ A U 



(49) 
(50) 

(51) 

(52) 



because of weak dependency of phase shift to the cosmo- 
logical model usually another model independent parameter 
which is so-called shift parameter 1Z from CMB observation 



11 oc 



iflat 

n 

ll ' 



(53) 



are used as another observational test. Where l[ ~" corre- 



j^flat 

sponds to the flat pure-CDM model with Q m = 1.0 and the 
same uo m and o»6 as the original model. It is easily show n 
that shift parameter is as follows <| J. R. Bond el. alJll997h 



1Z : 



-D L (z de c, f^roi ^A, p) 
(1 + Zdec) 



(54) 



The observational results of CMB experiments correspond 
to a shift parameter of 1Z = 1.716 ± 0.062 (given 
by WMAP, CBI, ACBAR) dD. N. Spergel el. al] l2003al : 
IT. J. Pearson el. all 120031 : IC L. Kuo et. al.ll2004l ). One of 



the advantages of using the parameter 1Z is its independency 
of Hubble constant. In order to put constraint on the model 
from CMB, we compare the observed sh ift parameter with 
that of model using likelihood statistic as ll J. R. Bond el. al.l 
1 19971 : lA. Melchiorri el. alJliool : I C. J. Odman et. al. Il2003h 



*cmb/ 2 



where 



2 _ [R-ob: 

XCMB — 



2 

°CMB 



(55) 



(56) 
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where 1Z t h and 1Z h s are determined using equation (|54|) 
and given by observation, respectively. Figure [10] shows con- 
stant value of I a in the joint space parameters (fl m , Qa) and 
(fi m , Q,x) for the braneworld and the ACDM model, respec- 
tively. Increasing (decreasing) Qa (&x) leads to an increas- 
ing in the value of present matter density to make constant 
value for I a- What we found is in agreement with figure [2] 
Another robust observational approach to inves- 
tigate cosmological models is inferring the behavior 
of the matter power spectrum and time evolution of 
gravitational clustering in both linear and nonlinear 
regimes. The simplest things to do are solving the rel- 
evant Boltzmann and Einestian eq uations for variou s 
matter contents in the Universe (|S. Dodelsonl [2003). 
Matter power spectrum and other non-linear effects can 
be a special tools to discriminate various models as well 
as to make more confined acceptable r; 



free parameters (see 


(IC. P. Ma et. all Il999l: 


Z. Mai 120061: 


T. koivistol2006l; G. Olivares et. all 2005, 2006; 


S. lee et. al. 


2006 


; V. R. Eke et. all 200ll; D. Jeong and E. Komatsu 


2006 


; U. Scliak and M. Zaldarriaea 


1999; 


D. H. Rudd et. al.ll2007l; A. Loeb and M. Zaldarriaea 


2005; 


R. R. R. Reis et. al. 


20051; G. Olivares et. al. 20051 


2006; 


P. J. E. Peebles! 2003 


: 1. P. Chimento and D. Pavonl 


2006; 


L. Amendolall2000l;lL. Amendola et. al.ll200ll; L. Amendola 


20021; L. Amendola et. al. 2003bl) for recent reviews). The 



convent ional form of mat ter power spectrum at the late 
time is |S. Dodelsonll2003h : 



u 

P{k,a) = 2ir 2 8 2 H T{kf 



Difl) 



(57) 



where n is the spectral index of the primordial adiabatic 
density perturbations, T(k) is transfer function determines 
the evolution of potential in the radiation-matter equality 
epoch and in the late time matter density fluctuations gov- 
ern by so-called growth function, D(a). 8h is also given by 
initial condition in the context of inflation, k is the wave- 
number of fluctuations in the Fourier space. Generally, it is 
well known that if any interaction between matter and dark 
energy in addition to the new kind of matter which are the 
responsible for the background dynamics of the Universe to 
be existed, they alter the matter power spectrum through 
the three following effects: The first one is that, the Hub- 
ble parameter in different models causes various dynamics 
for the background evolution (e.g. in our model is given by 
equation Q22[l ) as well as power spectrum. The second ef- 
fect is due to the inverse proportional of power spectrum 
to the matter density for a fixed potential, so any variation 
in the present value of matter density causes the smaller or 
larger amplitude for power spectrum. Third effect is related 
to the fact that in different cosmological models, the matter- 
radiation equality epoch, a eq and subsequently the value of 
fc e? change, so the turning over point in the power spectrum 
would be reformed. 

Here instead of observational constraint using mat- 
ter power spectrum we used the weakly model indepen- 
dent constraint by Baryon acoustic oscillation an d ig- 
nore any non-linear effects (|E. V. linden 120051 , [2003). Re- 
cently using the observations of large sc ale structures from 
the Slo an Digital Sky Survey (SPSS) lM~Tegmark, et~ai] 
11200481 ); |M. Tegmark. R. Michael BlantorTet all l|2004bl ) 
and Two Degree Field Galaxy Redshift Survey (2dFGRS) 



|M. Tegmark. et all (|2002l ). one can explore the validity of 
cosmological models. 

The large scale correlation function measured from 
46, 748 Luminous Red Galaxies (1RG) spectroscopic sam- 
ple of the SPSS includes a cle ar peak at about 100 Mpc 
IIP. J. Eiscnstci n el. al I 20051 ) . This peak was identified with 
the expanding spherical wave of baryonic perturbations orig- 
inating from acoustic oscillations at recombination. The co- 
moving scale of this shell at recombination is about 150Mpc 
in radius. In other words, this peak has an excellent match 
to the predicted shape and the location of the imprint of the 
recombination-epo ch acoustic oscillation on the l ow-redshift 
clustering matter (|D. J. Eisenstein el. al.ir20 05). Recently 
E. V. Under has shown in detail som e systematic uncer- 
taint ies for baryon acoustic oscillation (|E. V. linderl 120051 . 
2003). Nonlinear mode coupling related to this fact that 
ever though baryon acoustic oscillation is mostly contributed 
by linear scale, but the influence of non-linear collapsing 
has quite broad kernel. In other words, one might say that 
baryon acoustic oscillation are 90 — 99% linear in comparison 
to the CMB which is 99.99% linear, so this difference may 
affect on various models in different way. Careful works to 
constrain on the free parameters of underlying model needs 
to be carried out the effect of non-linear mode coupling in 
the results of constraint by SPSS observation. Nevertheless, 
roughly speaking regards the acceptance intervals for free 
parameter cover the real intervals determi ned by assum- 
ing nonlinearity mode for SPSS observation (| W. Hu et. al 



P. J. Eiscnstcin et. al. 1998; M. Amarzguioui et. al 
P. Eisenstein et. al. Il2004h . 



A dimensionless and independent of Ho version of SPSS 
observational parameter is 



A = D v (z sdss ) 



~sds 



H(z s d SB ;0. m , £Ia, Li)z 2 dss (l + z a d ss ) 



1/3 



(58) 

I 



where Dv(zsdss) is characterist ic distance scale of the surve. 
with the mean of redshift z a d«« |P. J. Eisenstein el. al.ll2005. . 
IC. Blake el. "aill2003l ; [S~Nesseris el. al.ll2007l ). We use the ro- 
bust constraint on the braneworld model using the value of 
A = 0.469 ±0.017 from the 1RG observation at z BdBB = 0.3 5 
|P. J. Eisenstein el. all I 20051 ; IP. Eisenstein et. al. 1 12004| ). 
This observation permits the addition of one more term in 
the x 2 °f equations (|44[) and (|56[) to be minimized with re- 
spect to H(z) model parameters. This term is 



2 |-4 bs 
XSDSS — 2 



■At: 



(59) 



-'SDSS 



This is the third observational constraint for our analysis. 

In the rest of this subsection we perform a combined 
analysis of SNIa, CMB and SPSS to constrain the param- 
eters of the braneworld model by minimizing the combined 
X 2 = XsNia+XCMB+xlDSS- The best values of the model pa- 
rameters from the fitting with the corresponding error bars 
from the likelihood function marginalizing over the Hub- 
ble parameter in the multidimensional parameter space are: 



n m = 0.29t°'°l Qa = -0.71+2;^ and fi = 



-0.021 

lcr confidence level with 



0.40+°™ at 



+0.00t° ;°* 



0.03 
n/Nd.o.f 

The Hubble parameter corresponding to 



Xmir* 



0.26 

0.93 demonstrated 
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Figure 11. Joint confidence intervals of f2 m and Qa> fitted with 
SNIa new Gold sample+CMB+SDSS. Solid line, dashed line and 
long dashed line correspond to 3a, 2a and la level of confidence, 
respectively. 



Figure 13. Joint confidence intervals of fi and Q m , fitted with 
SNIa new Gold sample+CMB+SDSS. Solid line, dashed line and 
long dashed line correspond to 3a, 2a and la level of confidence, 
respectively. 
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Figure 12. Joint confidence intervals of p, and Qa, fitted with 
SNIa new Gold sample+CMB+SDSS. Solid line, dashed line and 
long dashed line correspond to 3a, 2a and la level of confidence, 
respectively. 
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Figure 14. Joint confidence intervals of f2 m and Qa, fitted with 
SNIa SNLS+CMB+SDSS. Solid line, dashed line and long dashed 
line correspond to 3a, 2a and la level of confidence, respectively. 



the minimum value of % 2 1S Ho = 62.72. Here we obtain 
an age of 14.82+q'44 Gyr for the Universe (see section V 
for more details). Using the SNLS data, the best fit values 
of model parameters are: Q m = 0.27±o'.o2 = -0-7^to.°t 
and fi = O.OOlg'gQ at la confidence level with xLn/ft.o./ = 
0.86, states Qk = — 0.01 ^q'qI . Age of Universe calculating 
with the best fit parameters is 14.05^0 45 (see next section). 
Tables 2 and 3 give the best fit values for the free parame- 
ters and age of Universe computing with these values. Joint 
confidence intervals in free parameter spaces are shown in 
figures EH! 



5 AGE OF UNIVERSE 

The age of Universe integrated from the big bang up to now 
in terms of free parameters of the braneworld model is given 



by 

/•to 

io(fi m , Q.a, a*) = / dt 
Jo 

Figure [17] shows the dependency of Hoto (Hubble parameter 
times the age of Universe) on Qa and /j, for a flat Universe. 
Obviously increasing and /i result in a shorter and longer 
age for the Universe, respectively. As a matter of fact, ac- 
cording to the equation (|22[1 . Q,a behaves as inverse role of 
dark energy in the ACDM scenario and \i has the inverse 
role of w in the ACDM (see Figures [17] and [T8j ■ 

The "age crisis" is one the main reasons of the ac- 
celeration phase of the Universe. The problem is that the 
Universe's age in the Cold Dark Matter (CDM) Universe 
is less than the age of old stars in it. Studies on the old 
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Figure 15. Joint confidence intervals of ii and Qa, fitted with 
SNIa SNLS+CMB+SDSS. Solid line, dashed line and long dashed 
line correspond to 3a, 2a and la level of confidence, respectively. 



0.4 - 



0.35 - 



0.15 



I I I I ] I J | I J 




-0.5 







0.5 



Figure 16. Joint confidence intervals of /i and O m , fitted with 
SNIa SNLS+CMB+SDSS. Solid line, dashed line and long dashed 
line correspond to 3a, 2a and la level of confidence, respectively. 



stars (|E. Carretta el. aLll2000h su ggest an age of 13t% Gy r 
for the Universe. Richer et. al. (IH. B. Richer el. all l2002h 
and Hasen et. al. ||B. M. S. Hansen el. al.l I2002T ) also pro- 
posed an age of 12.7 ± 0.7 Gyr, using the white dwarf cool- 
ing sequence method (for f ull review of the cosmic age see 
(|D. N. Spergel el. al.|[2003ah ). Table 3 shows that age of the 
Universe from the combined analysis of SNIa+CMB+SDSS 
is 14.82to.44 Gyr and 14.05±g g Gyr for new Gold sam- 
ple and SNL S data, respectively, whil e ACDM implies 
13.7±0.2Gyr (|P. N. Spergel el. al.ll2003al]. These values are 
in agreement with the age of old stars (|E . Carretta el. al.l 
l200d : lL. M. Krauss et. alj|200ll ; [rTchabover et. al. ||2002h . 

To do another consistency test, we compare the age 
of Universe derived from this model with the age of old 
stars and Old High Redshift Galaxies (OHRG) in various 
redshifts. Here we consider three OHRG for comparison 
with the braneworld model, namely the LBDS 53W091, 
a 3. 5 - Gyr old radio gala xy at z = 1.55 Dunlop el. al.l 
1 19961 ; IH. Spinrardl 1 19971 ). the LBDS 53W069 a 4.0-Gyr 
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_1 34+17-34 


+0.20 
-0.23 


-0.76 


+0.11 
-0.16 


— 26 +1 ' 72 
U - ZD -1.00 


+0.27 
-0.23 


-0.76 


+0.16 
-0.17 


-0.2^ 


+0.02 
-0.02 


-0.74 


+0.04 
-0.02 




+0.05 
-0.04 


-0.74 


+0.07 
-0.05 


0.00«; 6 6 « 



Table 2. The best fit values for the parameters of the model 
using SNIa from new Gold sample and SNLS data, SNIa+CMB 
and SNIa+CMB+SDSS experiments at one and two a confidence 
level. 



Observation 



Age (Gyr) 



SNIa(new Gold) -0.26 

SNIa(new Gold)+CMB +0.01 

SNIa(new Gold)+ +0.00 
CMB+SDSS 

SNIa (SNLS) -0.90; 

SNIa(SNLS)+CMB -0.0] 

SNIa(SNLS)+ -0.0J 
CMB+SDSS 



-0.33 
-1.44 



+0.09 
0.10 



-0.04 
-0.04 



-1.64 
-0.59 



+0.23 
0.28 



+0.04 
0.03 



+2.10 



13.44 



13.53±i^ 
14-82±°-j5 



14.381?;™ 



14.38 



-3.03 
-14.38 



14 05+ 43 



Table 3. The best values for the curvature of the brane model 
with the corresponding age for the Universe from fitting with 
SNIa from new Gold sample and SNLS data, SNIa+CMB and 
SNIa+CMB+SDSS experiments at one and two a confidence 
level. 



old radio galaxy at z — 1.43 Dunlopl Il999l ) and a 
quasar, APM 8279 + 5255 at z = 3.91 with an age of 
t = 2 .ll° ? Gyr (IG. Hasinger el. al.ll2002i ; IS. Komossa et. all 
2002). The later has once again led to the "age crisis". An 
interesting point about this quasar is that it cannot be ac- 
commodated in the ACDM model (|D. Jain, el. al.ll2005l ). In 
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LBDS 


LBDS 


APM 


Observation 


53W069 


53W091 


08279 + 5255 




z = 1.43 


z = 1.55 


z = 3.91 


SNIa (new Gold) 


i nn+° 10 

±.uu_ 75 


1 0fi+ 010 
l.UO_ .7g 


65+ 14 
U - DO -0.41 


SNIa(new Gold)+CMB 


0.98+Hn 

—0.10 


— 0.10 


0.65+0-" 

— 0.07 


SNIa(new Gold)+CMB 


-, 99 +0.06 
1 - ZZ -0.04 


1 O1+0.06 

± -' 31 -0.05 


U ' 84 -0.04 


+SDSS 








SNIa (SNLS) 


9 qo + 1 ■ 15 

z.oo_ 1 13 


9 cc+1-35 
z - 00 -1.31 


1 O1+1-70 
- L - y± — 1.65 


SNIa(SNLS)+CMB 


i no+0.31 
l.zo_ 1 23 


, 09+O.34 
1 - oz -1.32 


n qk+0.28 
U - 85 -0.85 


SNIa(SNLS)+CMB 


1 16+ ' 04 


1 24+ - 04 
x - zyi -0.05 


79+°- 16 
u - ' 3 -0.04 


+SDSS 










Table 4. The value of r for three high redshift objects, using 
the parameters of the model derived from fitting with the obser- 
vations. 



order to quantify the age-consistency test we introduce the 
expression r as: 



tab 



tobsHo 



(61) 



where t(z) is the age of Universe, obtained from the equation 
(|17p and t b s is an estimation for the age of old cosmological 
object. In order to have a compatible age for the Universe we 
should have r > 1. Table 4 reports the value of r for three 
mentioned OHRG with various observations. We see that the 
parameters of braneworld model from the combined observa- 
tions provide a compatible age for the Universe, compared to 
the age of old objects, also in addition, SNLS data result in a 
shorter age for the Universe. Once again for the braneworld 
model, APM 08279 + 5255 at z = 3.91 has a longer age than 
the Universe but gives better r esult than most cosmologi- 
cal models investigated before (I L. F. Miranda et al.l l200ll ; 
IS. Rahvar elTal]|2007t Id. Jain, el. al.ll2005l ). 



0.95 




0.85 



Figure 17. -Hoio (age of Universe times the Hubble constant 
at the present time) as a function of Qa (upper panel) for a flat 
Universe and typical value of fi = —0.40. Increasing Qa gives a 
shorter age for the Universe. Lower panel shows the same function 
versus fi for the case Q m = 0.30, Qa = —0.70 (flat Universe). 



6 CONCLUSIONS AND DISCUSSIONS 

The impressive amount of data indicating a spatially flat 
Universe in accelerated expansion has posed the problem 
of dark energy and stimulated the search for cosmological 
models which are able to explain such unexpected behav- 
ior. Many rival theories have been proposed to solve the 
puzzle of the nature of dark energy ranging from a rolling 
scalar field to a unified picture where a single exotic fluid 
accounts for the whole dark sector (dark matter and dark 
energy). Moreover, modifications of the gravity Lagrangian 
have also been advocated. Although deeply different in their 
underlying physics, all these scenarios share the common fea- 
ture of well reproducing the available astrophysical data. On 
the other hand, alternative cosmology from the braneworld 
models provide a possible mechanism for the present accel- 
eration of the Universe congruously suggested by various 
cosmological observations. 

In braneworld scenarios, due to the usual energy con- 
servation law on the brane, we do not have energy flow from 



the brane onto the bulk or vice versa. There are numerous ef- 
forts to constrain the braneworld models but in all of them, 
there is no energy exchange between the brane and bulk. 
Theoretically, there are no fundamental reasons to forbid 
the energy exchange between the brane and bulk in a brane 
scenario. One can get this profile by relaxing the conserva- 
tion law on the brane. This energy exchange can alter the 
profile of the cosmic expansion and leads to a behavior that 
would resemble the dark energy. In this paper we focused 
our attention on the RS II braneworld model with energy 
exchange between the brane and bulk. We got the modified 
Friedmann equation (|22[) on the brane which can explain 
the cosmological behavior and describe a physically origin 
for the dark energy which is in good agreement with obser- 
vations. We explored the consistency of this scenario with 
the implication of up-to-date luminosity of supernova type 
la observed by two independent groups, new Gold sample 
and SNLS data set, acoustic peak in the cosmic microwave 
background anisotropy power spectrum and baryon acoustic 



14 M. Sadegh Movahed and Ahmad Sheykhi 





w 



Figure 18. Hoto (age of Universe times the Hubble constant at 
the present time) versus Q\ (upper panel) in the ACDM model 
for the case Qk = 0.0 and w = —1.0. Lower panel shows Hoto as 
a function of present equation of state, w, in the ACDM model 
with the typical values Cl m = 0.30 and Q\ = 0.70. 



oscillation measured by Sloan Digital Sky Survey. The effect 
of model free parameters on the matter power spectrum and 
the exploration of matter and dark energy interaction will 
investigate in our forthcoming paper. 

The best parameters obtained from the fitting with the 
new Gold sample data combined with CMB and SDSS ob- 



servations are: Q m = 0.29+o'„2, Qa = 



-0.7V 



and jx = 



— 0. 401^26 at l* 7 confidence level with Xmin/Nd.o.f = 0.93 
expressing spatially fiat Universe with Q.k = +0.00t 'o4- 



SNLS SNIa+CMB+SDSS give: 0, 



0.27 



—0.71 

Xmin 



+ 0.04 



and /i 



0.00 



+0.30 



+0.02 
■0.02 



at lcr confidence level with 

,+0.04 



/Nd.o.f = 0.86, asserting flu 
known ACDM mod e l impl ying -0.06 < Q, K < +0.02 
jD. N. Spergel el. all l2003al ) and some other interesting 
models such as Dvali-Gabadadze-Porrati (DGP) which in- 

using Gold 



O.Olt^ and fl K 



0.01 



+0.04 



dicates Ojr 

sample and SNLS data, respectively (IZong-Kuan Guo el. al.l 
l2006l : iMovahed et~al]|2007t iMovahed and Gassemill2007M . 

We also performed the age test, comparing the age of 



old stars and old high redshift galaxies with the age derived 
from this model. From the best fit parameters of the model 
using new Gold sample and SNLS SNIa, we obtained an age 
of 14.82to.44 Gyr and U.Q5±oaI Gyr, for the Universe, re- 
spectively. These results are in agreement with the age of 
the old stars. The age of Universe in this model is larger 
than what given in the other mod e ls (ID. N. Spergel el. aU 



L. F. Miranda et al.l l200ll : Is. Rahvar el. all 12007 
Movahed et. al.ll2007f ) 



2003al: 



To check the age crisis in this model we chose two high 
redshift radio galaxies at z = 1.55 and 2 = 1.43 with a 
quasar at z = 3.91. Two first objects were consistent with 
the age of Universe, i.e., they were younger than the Uni- 
verse while the third one was not but our model gave the 
better result than ACDM and a class of Quintessen ce model 
l| L. F. Miranda et al.|[200ll : IS. Rahvar el. al.ll2007l) . 

Finally, it must point out that the energy exchange term 
Q.A plays a crucial role in our work. In other words, in the RS 
II model without energy exchange where we have Qa = 0, 
we can not get late time acceleration expansion profile for 
our Universe! So we conclude that the usual RS II model 
should ruled out from present observational data. 
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